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Abstract 

In this paper, the non-Markovian dissipative dynamics of the phase 
damped Jaynes-Cummings model in the presence of a classical homo- 
geneous gravitational field will be analyzed. The model consists of 
a moving two-level atom simultaneously exposed to the gravitational 
field and a single-mode traveling radiation field in the presence of a 
non-Markovian phase damping mechanism. First, the non-Markovian 
master equation for the reduced density operator of the system in terms 
of a Hamiltonian describing the atom-field interaction in the presence of 
a homogeneous gravitational field will be presented. Then, the super- 
operator technique will be generalized and an exact solution of the 
non-Markovian master equation will be obtained. Assuming that ini- 
tially the radiation field is prepared in a Glauber coherent state and the 
two-level atom is in the excited state, the non-Markovian effects on the 
temporal evolution of collapses and revivals of the atomic population 
inversion and photon counting statistics of the radiation field in the 
presence of both the phase damping and a homogeneous gravitational 
field will be investigated. 
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1 Introduction 



Non-Markovian effects have received special attention in the past years, 
mainly in optics and radiation-matter interaction subjects, either in pre- 
dicting novel effects or due to the necessity to go beyond the Markovian 
approximation in experiments involving femtosecond processes. Among the 
experimental papers the recent ones can be cited. Tchenio et al. prepared 
a Non-Markovian atomic excitation process, with adjustable memory time, 
using correlated laser pulses and they verified that under strong-field condi- 
tions the atoms are not able to keep memory of the field phase and amplitude 
over a time interval larger than the coherence time [1] . Considering femtosec- 
ond experiments, non-Markovian behavior appears in the optical dephasing 
of molecules in solution, since the dynamics of the thermalized environment 
may occur on the same time scale of the system [2-6]. Concerning the theo- 
retical approach, Lewenstein et al. predicted the suppression of spontaneous 
emission related to the decay of cavity atoms in the presence of a strong 
driving field, thus modifying the spectrum of resonance fluorescence [7]. 
Villaeys et al. studied the non-Markovian effects in the atomic absorption 
band shape for the transient and steady-state regimes; they conclude that 
in the steady state the appearance of the non-Markovian effects are washed 
out and therefore they cannot be probed, but in the transient regime these 
effects are perceptible [8]. In this same line Gangopadhyay and Pay con- 
structed a non-Markovian master equation by considering density matrices 
with small delay time [9]. The prototype of such systems, proposed by the 
Jaynes-Cummings model (JCM) in 1963, [10] describes a two-level atom res- 
onantly interacting with a single-mode quantized field. It has proved to be 
a theoretical laboratory of great relevance to many topics in atomic physics 
and quantum optics, as well as in the ion traps, cavity QED and quantum 
information processing [11]. When the rotating wave approximation (RWA) 
is made, the model becomes exactly solvable and its dynamical features can 
be analytically brought to light revealing remarkable properties [12]. In the 
standard JCM, the interaction between a constant electric field and a sta- 
tionary (motionless) two-level atom is considered. With the development in 
the technologies of laser cooling and atom trapping the interaction between 
a moving atom and the field has attracted much attention [13-22] . 

Experimentally, atomic beams with very low velocities are generated in 
laser cooling and atomic interferometry [23]. It is obvious that for atoms 
moving with a velocity of a few millimeters or centimeters per second for a 
time period of several milliseconds or more, the influence of Earth's accel- 
eration becomes important and cannot be neglected [24]. A semi-classical 
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description of a two-level atom interacting with a running laser wave in 
a gravitational field has been studied [25,26]. However, the semi-classical 
treatment does not permit us to study the pure quantum effects occurring in 
the course of atom-radiation interaction. Recently, within a quantum treat- 
ment of the internal and external dynamics of the atom, we have presented 
[27] a theoretical scheme based on an su(2) dynamical algebraic structure to 
investigate the influence of a classical homogeneous gravitational field on the 
quantum non-demolition measurement of atomic momentum in the disper- 
sive JCM. Also, the effects of the gravitational field on quantum statistical 
properties of the lossless [28] as well as the phase-damped JCMs [29] were 
investigated. We reach to the point that the gravitational field seriously 
suppresses non-classical properties of both the cavity-field and the moving 
atom. Also, the effects of the gravitational field on the dynamical evolution 
of the cavity-field entropy and the creation of the Schrodinger-cat state in 
the Jaynes-Cummings model [30] are examined. 

On the other hand, over the last two decades much attention has been 
focused on the properties of the dissipative variants of the JCM. The theo- 
retical efforts have been stimulated by experimental progress in cavity QED. 
Besides the experimental drive, there also exists a theoretical motivation to 
include relevant damping mechanism to JCM because its dynamics becomes 
more interesting. A number of authors have treated the JCM with dissipa- 
tion by the use of analytic approximations [31,32] and numerical calculations 
[33-37]. The solution in the presence of dissipation is not only of theoretical 
interest, but also important from a practical point of view since dissipation 
would be always present in any experimental realization of the model. How- 
ever, the dissipation treated in the above studies is modeled by coupling 
to an external reservoir including energy dissipation. As is well known, in 
a dissipative quantum system, the system loses energy by creating a bath 
quantum. In this kind of damping the interaction Hamiltonian between bath 
and system does not commute with the system Hamiltonian. In general this 
leads to a thermalization of the system with a certain time constant. There 
are, however other kinds of environmental coupling to the system, which do 
not involve energy exchange. In the so-called phase damping [38] the inter- 
action Hamiltonian commutes with that of system and in the dynamics only 
the phase of system state is changed in the course of interaction. Similar to 
standard energy damping the off-diagonal elements of the density matrix in 
energy basis decay at a given rate. The phase damping can well describe 
some unaccounted decay of coherences in a single-mode micromaser [39]. It 
has also been shown that phase damping seriously reduces the fidelity of the 
received qubit in quantum computers due to the induced decoherence [40]. 
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The phase damping in the JCM with one quantized field mode has been 
studied [41]. The influence of phase damping on non-classical properties of 
the multi-quanta two-mode JCM has also been studied [42]. It has been 
found that the phase damping suppresses non-classical effects of the cavity 
field in the JCM. However, all of the foregoing studies have been done only 
under the condition that the influence of the gravitational field is not taken 
into account. 

In this paper, the non-Markovian dissipative dynamics of the phase 
damped Jaynes-Cummings model in the presence of a classical homogeneous 
gravitational field will be analyzed. The model consists of a moving two- 
Icvcl atom simultaneously exposed to the gravitational field and a single- 
mode traveling radiation field in the presence of a non-Markovian phase 
damping mechanism. In sect. 2, the non-Markovian master equation for the 
reduced density operator of the system in terms of a Hamiltonian describ- 
ing the atom-radiation interaction in the presence of a gravitational field 
will be presented. This Hamiltonian has been obtained based on an su(2) 
dynamical algebraic structure in the interaction picture. In sect. 3 an exact 
solution of the JCM with the phase damping in the presence of a gravita- 
tional field will be obtained, by which the dynamical evolution of the system 
is investigated. In sect. 4 the non-Markovian effects on both the cavity-field 
and the atomic properties in the presence of both the phase damping and a 
homogeneous gravitational field will be studied. Considering the field to be 
initially in a coherent state and the two-level atom in the excited state, the 
temporal evolution of the atomic inversion and photon counting statistics 
will be explored. Finally, our conclusions will be summarized in section 5. 



2 Non-Markovian Master Equation for the Phase 
Ddamped JCM in the Presence of a Gravita- 
tional Field 

The equation of motion for the density operator of the atom-radiation system 
and reservoir, Psr{t), in the Schrodinger picture is given by [29] 



dt 



-i[HT,Psr{t)]{n=l), 



(1) 



where 



HT = Hs + Hr + Vs 



(2) 
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with the Hamiltonian of the reservoir 



Hr = J2u;iblbi, (3) 

i 

and with the Hamiltonian of the interaction between the system and reser- 
voir 

3 

Vsr = HsY.Fj, (4) 



where 



6i and 6^ are the boson annihilation and creation operators for the reservoir 
and Ki is the coupling constant. The Hamiltonian Hg in (2) for the atom- 
radiation system in the presence of a classical gravity field with the atomic 
motion along the position vector x and in the RWA is given by {h = 1) 

X[exp{—iq.x)a) a- + exp[iq.x)a+a\, (6) 

where a and denote, respectively, the annihilation and creation operators 
of a single-mode traveling wave with frequency Uc, 9* is the wave vector of 
the running wave and a± denote the raising and lowering operators of the 
two- level atom with electronic levels |e), \g) and Bohr transition frequency 
LOeg- The atom-field coupling is given by the parameter A and p, x denote, 
respectively, the momentum and position operators of the atomic center of 
mass motion and g is Earth's gravitational acceleration. It has been shown 
[29] that based on an su(2) algebraic structure, as the dynamical symmetry 
group of the model, in the interaction picture, Hamiltonian (6) takes the 
following form 

hI= u;,{a^a + ^) + ^A$,g,t)So (7) 



+ 



where the operators 

-50 = ^(|e)(e| - \g){g\),S+=a\e){g\^,S. = ^b)(e|at, (8) 

2 Vk Vk 
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with the following commutation relations 

[So, S±] = ±S±, S\] = -2So, (9) 

are the generators of the su(2) algebra, the operator K = a^a + |e)(e| is a 
constant of motion which represents the total number of excitations of the 
atom-radiation, k{t) is an effective coupling coefficient 

^(t)=Aexp(-(A(#,i,5) + ^)), (10) 

and the operator 

MP:'tJ)=^c-{(^eg + J^+q-9t+—), (11) 

has been introduced as the Doppler shift detuning at time t [29]. By using 
following the same procedure as in refs. [29,44] we obtain the non-Markovian 
master equation for the reduced density operator of the system with neglect- 
ing 27ri^[J(c(;)|K(a;)p]|t^=o and the lamb shift term 

^ = -i[Hl, psit)] - ^[Hl, [Hi, mW (12) 

where Hg is given by (7). In Eq.(12), 7 and rj are the damping and the 
non-Markovian parameters, respectively, which depend on the temperature 

7 = 2.T/zm.-.o( ^^"^'"^"^l\ (13) 

OJ 

and 



77 = 2iTTp 



Jo ou 



where J{uj) and k{uj) are the spectral density of the reservoir and the cou- 
pling coefficient, respectively and p is the Cauchy principal part of the in- 
tegration [44]. 



3 Dynamical Evolution of the Non-Markovian Phase 
Damped JCM in the Presence of Classical Grav- 
ity 

In section 2, we reached to the non-Markovian master equation for the re- 
duced density operator of the atom-radiation system in the presence of a 



6 



classical homogeneous gravitational field. In this section, we now start to 
find the exact solution for the density operator Ps{t) of the non-Markovian 
master equation (12) with the Hamiltonian (7). For this purpose, the ap- 
proach presented in refs. [43-45] is applied. The formal solution is given by 

Ps{t) = exp(i?t) exp(S'i) exp{ft)ps{0), (15) 

e^p{Rt) = exp(i?it) exp(i?2i) exp(i?3t), (16) 

exp(f t) = exp(fit) exp(f2t), (17) 

where Ps(0) is the density operator of the initial atom-field system. The 
auxiliary super-operators Ri,R2,R3, S and Ti,Ti are defined through their 
action on the density operator such that 

eMRit)m = f: ^(^>p,(0)(l>, (18) 

eMR2t)psiO) ^ E ^-^^{H^fpsmHsf, (19) 
eMR3t)Ps{0) = E ^—^{H.rpMiH.f"', (20) 

ex.p{St)ps{0) = exjp{-iH^t)ps{0) e^piiHj), (21) 
exp(fit)p,(0) = eM-l{Hsft)ps{0) exp(-7(^')2t). (22) 

eMT2t)Ps{0) = exp{--f7]{H,)H)ps{0) cxp{--f7]{H,)H). (23) 

It is assumed that initially the radiation field is in a coherent superpo- 
sition of the Fock states, the atom is in the excited state |e), and the state 
vector for the center-of-mass degree of freedom is |V'c.m(0)) = / d^p(f)(j))\p}. 
Therefore, the initial density operator of the atom-radiation system reads 
as 

Ps{0) = Pfield{^) ® Patom{0) Pc.m(O) = 

where 

PfieldiO) = J212M0)wmmn){m\, (25) 

n m 

Pernio) =Jd''pJ dyr{p')<P{p)\p){p'\, (26) 



PfieldiO) (8) Pernio) 





(24) 
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with Wn{0) 



cxp(- 



-. The Hamiltonian (7) can be expressed as a sum 
of two terms which commute with each other, that is, 



H,=Hi + H2,[Hi,H2]=0 



where 



H2 = -A{p,g,t)So + {k{t)\/KS- + k*{t)\/KS+). 
In the two-dimensional atomic basis we have 



Hi = iOc 



n + 



n 



Ho 



(27) 

(28) 
(29) 

(30) 
(31) 



K(t)aJ 4 

Also, the square of the Hamiltonian (7) can be expressed as a sum of two 
operators, one of them is diagonal, in the form 

I 



where 



Hf + Hl 



a;2(n + i)2 + A2(n + l) + (^^) 



Ai + A2, 



(32) 



(33) 



a;2(n-i)2 + A2n+(^^^)2 



and 



A2 = 2H1H2 



{n + U^) {n+\)n*{t)a 
{n-\)n{t)d) -(n-i)(^ 



(34) 



It is easily proved that [741,^2] = 0. Similarly, the square of the (i^,)^ can 
be expressed as a sum of two operators, one of them is diagonal, in the form 



{H,f = ^3 + ^4, 



where 



A^ = Al + Al 



(^3)11 

{A^)22 



(35) 
(36) 
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with 



(^3)11 



1. 



2i2 



+ 



[co-i{h + -f + x\h + i) + c- 4 

Au^n + ^)(^^) + 4a;^^(n + '-fin + 1), 



(37) 



(A 



3)22 



^-2 , x2^ , ,^(P>5,i).|2]2 



['^c(n-2)' + A^n + ( ^ 



(38) 



1. 
2^ 



and 



with 



Ai = 2A1A2 



(^4)11 (44)12 
(^4)21 (^4)22 



(39) 



(i4)n = [-e^(n+^)^ + A2(n+l) + (^^)^]4a;.(n+i)(^^Mll)), (40) 
(^4)12 = k'(n + ^)^ + \\h + 1) + (^^MlI) )2]4^^(^ + 1 (41) 
(i4)2i = [ul{n - \f + X'h + (^iM^)2]4u;,(n - \)Kal (42) 



(A4) 



22 



r 2/-. 1x2 , ^2^ , AiP^9,t) 2u I A{p,g,t) 
-['^c("-o) +An+( )]4u;c(n--)( ). (43) 



4 ^ ' 2'^ 4 

It is easily proved that [A3, A4] = 0. Taking into account the initial condition 
(24) the auxiliary density operator P2{t) is defined as 

P2(t) =eMSt)c^Y>{ft)ps{Q) (44) 
= cxp(—iH2t) cxp{—jA2t) exp{—jriA^t) 
xpi (t) exp(-7r/A4t) exp(-7A2t) exp{iH2t), 

where the operator pi (t) is defined by 

pi{t) = \^it)){^{t)\^\e){e\, (45) 

with 

Mt)) = exp(-7(l + m^'cin + If + \\n + 1) + {^^^f]) 



(46) 



xexp(-47??a;^(n + -)'(A'(n + l) + ( ^ 



Prom (30) and (33) we have, respectively 

exjp{—iujct{n + i)) 



exp{—iHit) ■■ 

ex.p(—^Ait) 









eyip{—iu;ct{h — i)) 



(^i)ii(n,t) 

(ii)22(n,i) 



(47) 

(48) 



where 



(ii)n(n, t) = exp(-7i[a;2(n + + A'(n + 1) + i^^^^f]), (49) 



iAi)22in, t) = eM-lt[^l{n - -f + A^n + ( ^ 



Also, we can write 

exp(-7^2i) 



exp(— 777^4^) 



ei{n,t) e2{n,t)a 
63(71, t)a) 64(71, t) 

e'i{h,t) e'2{n,t)a 



?])■ (50) 

(51) 
(52) 



where 



ei(ra,t) = cosh(7i^ci(n,t))-a;c( — ){n + -) 



MP,9,t).,^ , 1 sinh(7tA/ci(n,t)) 



2 ^/cl{n,t) 



,N o x/-^ 1 smh 7t^ci n- l,t)) 
e2(n,t) = -2uJcX{n - -) , 

2 v^ci(n-l,t) 



.X o 1. smh(7t^C2(n,t)) 

2 VC2{n,t) 



^{P^9,t)^,^ 1. sinh(7t^C2(n,i)) 



\/c2{n,t) 



, (53) 

(54) 
(55) 
, (56) 



with 



ci(n, t) = a;,^(^^)^(n + \f + X\^^^f{n + l)(n + i)^ (57) 



C2(n,t) = 



2AiP,9,t),2^, 1,2 , ,2,A(p,5,t)^2^^._1^2^ ^^g^ 



2 + 2 
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and 



'i(n,t) = cosh(7??tyc'i(n,t))-2a;c( ^ )("+2)^i("' *) > 

yc'i(fi,t) 

(59) 



(60) 



.X o xr-^ 1.^ , sinh(7r?tVc'i(ri-l,i)) 
e'2{n,t) = -2uJcX{n- -)L2{n,t)- 

^ 'c'i(n-l,0 



1 ^ sinh(7r7tvc'2(n,t)) 
e's{h,t) = -2u;eA(n- -)L2(n,t) ^7 ^ ^ (61) 



e'4 



r-^ +^ v,^ . /7 .A(^M^^' /-^ ,, sinh(7r?t7c^2(n,t)) 

(n,t) = cosh(7r?tyc'2(n,t))-a;c( 2 )('^-2^^^"' ^ > 

yc'2(n,t) 

(62) 

with 



where 



(63) 

2,A(p,5,t)^2(^_l)2^2(^^^) + ^^(^^^)^^(^^^)^(^_l)2^ (64) 



Li(n, t) = u^in + + A^(^ + 1) + (^^)^ 



2/-^ 1n2 , \2^ , AiP^9,t).2 



L2ih,t)=u;iin--y + X'h + { ^ 



(65) 



(66) 



Similarly, the operator exp{—iH2t) in the two-dimensional atomic basis can 
be stated as 

di{n,t) d2{n,t)a 
d^{n,t)a) di{n,t) 



exp{—iH2t) 



(67) 



where 

di{n,t) = cos(t(( 



A(p,5,t),2, ^2,-^ ,A(p,5,i),sin(t((^MM)2 + A2(n + l))) 



-)^+A^(n+l)))-(- 



-)- 



c?2(n, t) = —iX 



sin(t(( 



Mp,9,t) \2 I \2 



^(A(gM)2 + A2(n + l) 



4 ^ ^(^(M*))2 + A2(n + 1)) 

(68) 

(69) 



+ A^(n + 1))) 
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dsin, t) 



(70) 



d4{'h,t) = cos(t(( 



^^Mlllf + A^n)) - ( 



Then, from (51) and (67), it follows that 
exp(-z^2*) ex.v{-lMt) ■■ 

where 



fi{n,t) f2in,t)a 



h{n,t) = ei{n,t)di{n,t) + e2{n,t)d2{n,t), 
f2{n,t) = e2in,t)di{h,t) + ei{h,t)d2ih,t), 
h{n,t) = e3{h,t)di{n,t) + ei{n,t)d3{h,t), 
h{n, t) = 64(71, t) J4(n, t) + e3(n, t)d'i{h, t). 
Also, from (52) and (60), we have 



exp(— z/f2i) exp(— 7^2*) exp(— 777A4t) = 
where 



ji(n,t) j2{n,t)d 
J3{h,t)d^ J\{n,t) 



(71) 

(72) 

(73) 
(74) 
(75) 
(76) 

(77) 



Ji{n,t) = fiih,t)e\ih,t) + f2in,t)e'2in,t), (78) 

J2{n,t) = f2{h,t)e'i{n,t) + fi{n,t)e'2{n,t), (79) 

Mn, t) = h{n, t)eU{n, t) + U{n, t)e'3{h, t), (80) 

Mh,t) = U{n,t)eU{n,t) + h{n,t)e'3{h,t). (81) 

Substituting (45) and (77) into (44), an explicit expression for the operator 
P2{t) can be obtained as follows 

(/52(^))^,,■ = |^^(^))(*,■(^)|,(^,i = l,2), (82) 

with 

|*i(t)) = Ji{n,t)\^{t)), |^2(t)) = Un,m{t)), (83) 

where |^'(t)) is given by Eq.(46). Now, we obtain the action of the operator 
exp(iii) = exp(i?ii) exp(i?2i) exp(ii3i) on the operator P2{t) 



i=0 



(84) 
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where 



fe=o 

k 



k\ 



which can be expUcitly expressed as follows 



Tjk — l Tjl 



where 



K{t) 



g+in,t) 

■0* {h,t) 



at 



^(^)2+A2(n+l)) 



gt{n,t) = vX{h,t)-^^P4^vt{h,t), 



with 



2 + A2(n + l)), 



s±{h,t) =Uc{n - -) ± W( 



A(p,5,t) 



)2 + A^n. 



(85) 
(86) 

(87) 



(89) 
(90) 



ui{h,t) = ^ifl{h,t)±ft{h,t)),vi{h,t) = ^{sX{h,t)±stih,t)), (91) 



(92) 



(93) 



Finally, by substituting (85) and (87) into (86) we obtain the exact solution 
of the non-Markovian master equation (12) for the phase damped JCM in 
the presence of a classical homogeneous gravity field 



Ps{t) 



where 



(/5.)ii(t) (/5.)i2(t) 

{Ps)2l{t) {Ps)22{t) 



k=0 



k\ 



(94) 



(95) 



13 



+ 
+ 
+ 



av'}inMk)22it)v}in,t)a^)\m\'', 



(27*) 



k 



{Ps)22{t) = E '-^{vJ^ih,t)a^{pA)u{t)av'J^{n,t) (96) 



A;=0 



A;! 



+ 9-{n,t){k)2i{t)avJ^{n,t) 

+ 9'-{nMM22m'i{n,mmf, 



iPsMt) = {Ps)2i{t)^ = E ^^{vJ^{n,t)a\p4)u{t)gl{h,t) (97) 

fe=o 



with 



where 



+ g'l{n,t){p4)2i{t)agl{h,t) 

+ v}{h,t)a^{Mi2{tmvj!{h,t)a^ 

+ 5-(n,t)(p4)22(t)^)l'(n,t)at)|</,(p)|2, 



v'}ih,t) = ■.v'_ih,t) (98) 

V(^^)2 + A2n 



(/54)ii(i) = E ^^^(5+'(n,i)(/$3)ii(i)5i(n,t) (99) 

+ avy{h,t){p3)2i{t)g'_^{h,t) 
+ i^{n,t){p^)i2{t)vl{n,t)a^ 
+ avy{h,t){p3)22{t)vl{h,t)a^), 

(P4)22(i) = E ^^^^(«-^(n,i)a^(/53)ii(0a^)^(^,*) (100) 
j=o 3- 

+ gj'{n,t){ps)2i{t)avl(fi,t) 

+ i;l2^'(n,t)at(p3)i2(t)£(n,t) 

+ £^(n,t)(p3)22(i)9'-(n,t)), 
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(P4)l2(i) 



+ 
+ 
+ 



(P4)2l(t)^ 

g- (n, t){p3)2i{t)ag]_{n, t) 

v'y (n, t)a^ (p3)i2 (t) v^' (n, t)at 

5-'(n,t)(/53)22(t)0^'(/i,i)a^), 



(101) 



with 



+ 
+ 
+ 



E Mr^ (5; * ) (P2 ) 11 i ) (102) 



i=0 



avi{n,t){p2)2i{t)g+{n,t) 



{h)22{t) = E ^r^(^>^(n,t)at(p2)ii(t)at)?(n,t) (103) 

1=0 



+ gUn,t){p2)2i{t)av^'{n,t) 
+ i;^(n,t)at(^2)i2(t)5-(n,t) 
+ gl{h,t)ip2)22{t)g^-^{n,t)), 

(/53)i2(i) = (/54)2i(t)^ (104) 

= E *)«^ (P2) 11 (n, i) 

+ 5L(n,t)(p2)2i(t)a5+(ri,t) 

+ v'l{h,t)aHp2)i2mt)v'^\n,t)a^ 

+ gl{h,t){p2)22{t)vy{h,t)a^), 

where we have defined {p2{t))i,ji {hi = 1,2) in (82). 

Making use of the solution given by (94), one can evaluate the mean values 
of operators of interest. In the next section, it will be used to investigate 
various dynamical properties of the non-Markovian phase damped JCM in 
the presence of a homogeneous gravitational field. 
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4 Dynamical Properties 



In this section, the non-Markovian effects on the quantum statistical prop- 
erties of the atom and the quantized radiation field in the presence of both 
the phase damping and the gravitational field will be studied. 



4a. Atomic Population Inversion 

The atomic population inversion is expressed by the expression 

Wit) = {asit)) = Tratomipatomit)^3it)), (105) 

where 

Patomit) =Trfield{Ps{t)- (106) 

It can be rewritten (105) as follows 

/« CO 
£p J2 {i\patom{t)a3{t)\i) = J d^pY,{{n\(^{{e\ps{t)\e)-{g\ps{t)\g))(^\n)). 
i=e,g n=0 

(107) 

Therefore, by using (94) and (107) we obtain 

W{t) = / clMT. E ^-^i{n\{Ps)uit)\n) - {n\{ps)22{t)\n))), (108) 

fe=0 n=0 

where {p2{t))i,j, (i, j = 1, 2) is given by (82). 

In Fig.l, the atomic population inversion as a function of the scaled 
time Xt for different values of the parameters q.g, 7 and i] are plotted. In 
this figure and all the subsequent figures we set q = Wm-^, M = W'^^Kg, 
g = 9.8^, Urec = llJ = -5 X 10"^^, A = 1 X 10^^:^, a = 2, A = 1.8 x lO^if^, 
~ s/i-KOQ ^-^P(~^) Co = 1 [25-30]. In Fig. la, three parameters: 
q.g = 0, 7 = and 77 = are considered. When q.g = 0, the momentum 
transfer from the laser beam to the atom is only slightly altered by the gravi- 
tational acceleration because the latter is very small or nearly perpendicular 
to the laser beam. When 7 = 0, there is no the phase damping and 77 = 
means that we consider Markovian approach. As it is seen from Fig. la for 
the atomic population inversion the Rabi-like oscillations can be identified. 
In Figs. lb and Ic, small gravitational influence q.g = 0.1 x 10'' with 77 = for 
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7 = and 7 = 7 x 10~^^ is considered, respectively. By considering both 
small gravitational influence and small non-Markovian effect in the presence 
of the phase damping the Rabi oscillations of the atomic population inver- 
sion disappear (see Fig. Id). By comparing Figs. Ic, Id and If, the influence 
of the non-Markovian on the time evolution of the atomic population inver- 
sion when there are both the phase damping and the gravitational field can 
be seen. 



4b. Photon Counting Statistics 



Now, the influence of gravity on the sub-Poissonian statistics of the ra- 
diation field will be investigated. For this purpose, we calculate the Mandel 
parameter defined by [46] 

For Q < {Q > 0), the statistics is sub-Poissonian (super- Poissonian); 
(5 = stands for Poissonian statistics. Since {n{t)) = Yl^=onP{n,t) and 
{n{tf) = En=on^P{n,t) we have 

00 00 00 

Qit) = {{[J2n'Pin,t)] - [Y,nPin,t)f}[J2nPin,t)]-') - 1, (110) 

n=0 n=0 n=0 

where the probability of finding n photons in the radiation field is found to 
be 

P{n,t) = {n\pfieid{t)\n) = {n\TratomPs{t)\n), (111) 

and by using (94) we have 

r 00 00 (r, .\k 
P{n,t) = d'pJ2Y. ^-^i{r^\iPs)ii{t)\n) + (n|(p,)22(t)|n)). (112) 

k=On=0 

Therefore, by using (110) and (112) we obtain 



00 . 00 00 (n f\k 

n=0 k=On=0 

x{{n\{ps)n{t)\n) + {n\{ps)22{t)\n)))] 



(113) 
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oo „ oo oo f\k 

Y.n{hfpY.Y. HA(("l(P3)ii(*)|n) + {nm22{t)\n)))f} 

re=0 fc=On=0 

oo „ oo oo (r^ f\k 

JLniJd'pY^Yl ^-^i{n\iPs)uit)\n) + (n|(p,)22W|n)))]-i) - 1. 



The numerical results for three values of the parameter q.g, 7 and rj 
are shown in Fig.2. As it is seen, the cavity-field exhibits alternately sub- 
Poissonian and super-Poissonian statistics when every three influences are 
negligible. For small gravitational influence and small non-Markovian ef- 
fect in the presence of the phase damping or with increasing rj the sub- 
Poissonian characteristic is suppressed and the cavity-field exhibits super- 
Poissonian statistics. After some time, the Mandel parameter Q is stabilized 
at an asymptotic zero value; the larger the parameter r] is more rapidly Q{t) 
reaches the asymptotic value zero. 



5 Summary and conclusions 

In this paper, the non-Markovian dissipative d3mamics of the phase damped 
Jaynes-Cummings model in the presence of a classical homogeneous gravi- 
tational field have been analyzed. The model consists of a moving two-level 
atom simultaneously exposed to the gravitational field and a single-mode 
traveling radiation field in the presence of a non-Markovian phase damp- 
ing mechanism. First, the non-Markovian master equation for the reduced 
density operator of the system in terms of a Hamiltonian describing the 
atom-field interaction in the presence of a homogeneous gravitational field 
has been presented. Then, the super-operator technique is generalized and 
an exact solution of the non-Markovian master equation is obtained. Assum- 
ing that initially the radiation field is prepared in a Glauber coherent state 
and the two-level atom is in the excited state, the non-Markovian effects on 
the temporal evolution of collapses and revivals of the atomic population 
inversion and photon counting statistics of the radiation field in the pres- 
ence of both the phase damping and a homogeneous gravitational field have 
been investigated. The results are summarized as follows: with increase of 
the non-Markovian parameter r] for small values of the damping parameter 
7 and gravity-dependent parameter q.g, 1) the Rabi-like oscillations in the 
atomic population inversion disappear and 2) the sub-Poissonian behaviour 
of the cavity-field is suppressed and it exhibits super-Poissonian statistics 
and after some time, the Mandel parameter Q{t) is stabilized at an asymp- 
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totic zero value; the larger the parameter 77 increases Q{t). 
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FIGURE CAPTIONS: 

FIG. 1 Time evolution of the atomic population inversion versus the 
scaled time At. Here we have set q = 10^ m~^, 

M = 10-26%,5 = 9.8f,Urec = -5 X W^'-f, 

A = 1 X 10^^, = 0, a = 2, A = 1.8 X 10^^, 

a) For q.g = 0, 7 = 0, = 0. 

b) For q.g = 0.1 x lO'^, 7 = 0, 77 = 0. 

c) For q.g = 0.1 X 10^ 7 = 7 X lO'^^, r? = 0. 

d) For q.g = 0.1 x 10^, 7 = 7 x 10'^^, ri = 5 x 10"^^. 
f)For ^.^= 0.1 X 10^, 7 = 7 X 10"^^, r? = 5 x 10"^^. 

FIG. 2 Time evolution of the atomic dipole squeezing versus the scaled 
time Xt with the same corresponding data used in Fig.l; 

a) For q.g = 0, 7 = 0, 77 = 0. 

b) For q.g = 0.1 x 10^, 7 = 0, r? = 0. 

c) For q.g = 0.1 X 10^, 7 = 7 X 10"^^, r] = 0. 

d) For q.g = 0.1 x lO'^, 7 = 7 x IQ-^^^f^, rj = 5 x 10'^^. 
f)For q.g = 0.1 x 10^, 7 = 7 x 10"^^, r) = 5 x IQ-^^. 



22 



